A common mode of origin of power laws in models of market and earthquake 
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We show that there is a common mode of origin for the power laws observed in two different 
models: (i) the Pareto law for the distribution of money among the agents with random saving 
propensities in an ideal gas-like market model and (ii) the Gutenberg-Richter law for the distribution 
of overlaps in a fractal-overlap model for earthquakes. We find that the power laws appear as the 
asymptotic forms of ever-widening log-normal distributions for the agents' money and the overlap 
magnitude respectively. The identification of the generic origin of the power laws helps in better 
understanding and in developing generalized views of phenomena in such diverse areas as economics 
and geophysics. 
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I. INTRODUCTION 

Recently we have shown that the Pareto law appears asymptotically (m — > oo) in the distribution of money m 
among the agents in the steady state of a trading market model: 

P{m) ~ m- {1+u) , v=l (1) 

when the agents have random saving propensities 0, Q • The market is modeled as an ideal gas where each molecule 
is identified with an agent, with the additional attribute that each agent has a random saving propensity, and each 
trading event between two agents is considered to be an elastic or money conserving collision between two molecules. In 
another model — a geometric model for earthquakes Q — we have shown that a power law similar to the Gutenberg- 
Richter law 0, appears in the asymptotic distribution of the overlap S between two dynamically intersecting Cantor 
sets: 

G(S) 7 = 1. (2) 

Since a geological fault is formed of a pair of fractal rock surfaces that are in contact and in relative motion, it is 
modeled by a pair of overlapping Cantor sets (the simplest known fractal), one shifting over the other. The overlap 
between the two Cantor sets represents the area of contact between the two surfaces of the fault and hence it is 
proportional to the energy released in an earthquake resulting from ruptures in the regions of contact. 

In both the models we get simple power laws with the exponents v = 7 = 1. Although these have been obtained 
separately for the two models, using both numerical and analytic methods, we show here that the two cases have a 
common feature that results in a common mode of origin of the power laws observed in the distribution of money 
m and fractal overlap S. The derivation of the power laws presented here shows that the common feature is a log- 
normal distribution in which the normal factor spreads indefinitely thus leaving the power-law factor to dominate the 
asymptotic distribution. 

II. THE IDEAL GAS MARKET MODEL 

Let us first consider the ideal gas model of an isolated economic system — that we refer to as the 'market' — in 
which the total money M and the total number of agents N are both constant; there is neither any production nor any 
destruction of money within the market and no migration of agents occurs between the market and its environment. 
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FIG. 1: Schematic diagram of the trading process. Agents i and j redistribute their money in the market: rrii(t) and rrij(t), 
their respective money before trading, changes over to rrii(t + 1) and rrijit + 1) after trading. 

The only economic activity allowed in the market is trading among the agents. Each agent i possesses an amount 
of money rrii(t) at time t. The time t is discrete and each event of trading is counted as a unit time step. In an 
event of trading (shown schematically in Fig. 1) a pair of agents i and j randomly redistribute their money between 
themselves such that their total money is conserved and none of the two agents emerges from the trading process 
with negative money (i.e., debt is not allowed): 

rrii(t + 1) + mj(t + 1) = m^t) + mj(t), (3) 
nii(t) > for all i at all t. 

It has already been shown that in the steady state market (t — * oo) the money m with the individual agents follow 
the Gibbs distribution Q: 

P( m ) = -exp(--J; T=— (4) 

when there is no restriction on the amount of money each agent can trade with except that it must satisfy the 
conditions of Eq. . Here T represents the economic equivalent of temperature and is defined as the average money 
per agent in the market. If each agent saves a fraction A (0 < A < 1) of its own money at every trading and A is the 
same for all agents at all time steps, the individual money with the agents in the steady state follows the Gamma 
distribution [8j. If we consider the effect of randomly distributed saving fraction A^ among the agents i, the money 
distribution in the steady state assumes the form of the Pareto law. The evolution of the agents' money in a trading 
can be written as 

m,(< + 1) = AjTOi(t) + e t [(1 - X l )m i {t) + (1 - Xj)mj(t)] (5) 

and 

ro,-(i + 1) = Xjrrijit) + (1 - e t ) [(1 - X^m^t) + (1 - Xj)mj(t)] (6) 

where Xi and Xj are the saving fractions of agents i and j respectively. The saving fractions Xi are quenched, i.e., 
fixed in time for each agent i and are distributed randomly and uniformly (like white noise) on the interval [0,1). 
The random division of the total traded money is given by the number < e t < 1 that varies randomly with the 
trading events t. The money distribution in the steady state is found to have a long power-law tail (shown in Fig. 2) 
that fits with the Pareto law for v = 1 Q . We also have analytic proofs 0, |tj 0] of the Pareto distribution of money 
observed in this random-saving gas- like model; all these proofs proceed by formulating the trading events as scattering 
processes and show that the Pareto distribution is a steady state solution of the scattering problem. 

Here we give a simple derivation of the asymptotic distribution of money in the steady state of the market model 
using an argument of the mean-field type, thereby avoiding the intricacies of the previous proofs. In our approach 
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FIG. 2: The money distribution P(m) for (a) random et and (b) et = 1/2. The power law in both cases has the same exponent: 
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Numerical simulation has been done for a system of N = 200 agents. 



the money redistribution equations JSJ and JBJ are reduced to a single stochastic map by taking the product of the 
two equations: 



m,i(t + l)rrij(t + 1) = aj(e t , Xi)m 2 (t) + otj(e t , Aj)m 2 (t) + <%-(e t , Ai, \j)mi(t)m,j(t). 



(7) 



Now we introduce a mean- field-like approximation by replacing each of the quadratic quantities to 2 , to 2 and mirrij 
by a mean quantity to 2 . Therefore Eq. J7J) is replaced by its mean-field- like approximation 



TO 2 (t+ 1) = r)(t)m 2 (t) 



(8) 



where ry(t) is an algebraic function of A^, Xj and e t ; it has been observed in numerical simulations of the model that 
the value of et, whether it is random or constant, has no effect on the steady state distribution 9] (illustrated in Fig. 
2) and the time dependence of rj(t) results from the different values of Aj and A^ encountered during the evolution of 
the market. Denoting log(m 2 ) by x, Eq. JSJ) can be written as: 



x(t + l) = x(t)+S{t), 



(9) 



where 8(t) — log?y(t) is a random number that changes with each time-step. The transformed map (Eq. |5J) depicts a 
random walk and therefore the 'displacements' x in the time interval [0, t] follows the normal distribution 



T{x) ~ exp ( — - 



Now 



V(x)dx = P(m)dm 2 



where P(m) is the log-normal distribution of m . 
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P(m) ~ — exp 
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(log(m 2 )Y 



(10) 



(11) 



(12) 



The normal distribution in Eq. (|10() spreads with time (since its width is proportional to y/i) and so does the normal 
factor in Eq. H12|l which eventually becomes a very weak function of m and may be assumed to be a constant as 
t — * oo. Consequently P(m) assumes the form of a simple power law: 



P(to) 
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oo, 



(13) 



that is clearly the Pareto law (|TJ for the model. 
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FIG. 3: The second generations (n = 2) of two Cantor sets with periodic boundary conditions, one shifting uniformly over the 
other. The total length of the shaded regions gives the overlap S between the two sets at any instant of time. The distribution 
of log 2 S is given by Eq. 1141 that is seen to approach a log-normal distribution of S (Eq. 1191 with increasing n and will 
eventually assume the asymptotic power law form (Eq. 1201 as n — > oo. 



III. THE FRACTAL-OVERLAP MODEL OF EARTHQUAKE 



Next we consider a geometric model A] of the fault dynamics occurring in overlapping tectonic plates that form the 
earth's lithosphere. A geological fault is created by a fracture in the earth's rock layers followed by a displacement 
of one part relative to the other. The two surfaces of the fault are known to be self-similar fractals. In this model 
a fault is represented by a pair of overlapping identical fractals and the fault dynamics arising out of the relative 
motion of the associated tectonic plates is represented by sliding one of the fractals over the other; the overlap S 
between the two fractals represents the energy released in an earthquake whereas log S represents the magnitude of 
the earthquake. In the simplest form of the model each of the two identical fractals is represented by a regular Cantor 
set of fractal dimension log 2/ log 3. This is the only exactly solvable model for earthquakes known so far. The exact 
analysis of this model for a finite generation n of the Cantor sets with periodic boundary conditions showed that 
the probability of the overlap S, which assumes the values S = 2 n ~ k (k = 0, . . . , n), follows the binomial distribution 
F of log 2 S = n — fc: 



Pr (S = 2 n - k ) = Pr (log 2 S = n-k) 

n — k / n\ k 

= F(n-k). 



n 

n — k 



(14) 



Since the index of the central term (i.e., the term for the most probable event) of the above distribution is n/3 + S, 
—2/3 < S < 1/3, for large values of n Eq. (|14fl may be written as 



F 



n 

n ± r 



2\^ 



(15) 



by replacing n—k with n/3±r. For r -C n, we can write the normal approximation to the above binomial distribution 
as 



F 



/ 27rn 



exp 



9r2 
' ~2n 



(16) 



Since log 2 S — n — k = § ± r, we have 



F(log 2 5)~ ^exp 



(log 2 S) 2 



(17) 
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not mentioning the factors that do not depend on S. Now 



F (log 2 S) d (log 2 S) = G(S)dS 



(18) 



where 



G(5)~ -exp 



(log 2 ^) 2 



(19) 



is the log-normal distribution of S. As the generation index n — > oo, the normal factor spreads indefinitely (since 
its width is proportional to y/n) and becomes a very weak function of S so that it may be considered to be almost 
constant; thus G(S) asymptotically assumes the form of a simple power law with an exponent that is independent of 
the fractal dimension of the overlapping Cantor sets: 



This is the Gutenberg-Richter law (Eq. [2J for the fractal-overlap model of earthquakes. It was also observed in 
numerical simulations [T^ | that G(S) ~ S 1-7 , 7 w 1 for several other regular and random fractals, thus suggesting 
that the exponent may be universal. 

The exact result of Eq. (|14l) in Ref. disagreed with the asymptotic power law of Eq. H2U|) obtained previously by 
renormalization group analysis of the model for n — ► 00 in Ref. [j] . The disparity between the two results had appeared 
because it was overlooked that the former is the exact distribution of log 2 S whereas the latter was the asymptotic 
distribution of S. However the above analysis shows that the power law in Eq. I|20|l is indeed the asymptotic form 
of the exact result. This is qualitatively similar to what is observed in the distribution of real earthquakes: the 
Gutenberg-Richter power law is found to describe the distribution of earthquakes of small and intermediate energies; 
however deviations from it are observed for the very small and the very large earthquakes. 

The fact that the fractal-overlap model produces an asymptotic power law distribution of overlaps suggests that 
the Gutenberg-Richter law owes its origin significantly to the fractal geometry of the faults. Furthermore, since this 
model contains the geometrical rudiments (i.e., the fractal overlap structure) of geological faults and it produces an 
asymptotic distribution of overlaps that has qualitative similarity with the Gutenberg-Richter law, we are inclined to 
believe that the entire distribution of real earthquake energies is log-normal that is wide enough for the Gutenberg- 
Richter power law to be observed over a large range of energy values. 



In the trading market model, we have shown that the money redistribution equations for the individual agents 
participating in a trading process can be reduced to a stochastic map in m 2 (Eq. (SJ. Using the transformation 
x = log(m 2 ), the map was reduced to a random walk in the variable x and hence the distributions of x and m 2 were 
found to be normal and log- normal respectively; in the steady state, i.e., for t — > 00, the latter was found to assume 
the form of a power law identical to the Pareto law with the exponent v = 1. Likewise, in the fractal-overlap model 
for earthquakes the distribution of overlaps was found to be log-normal for large generation indices n of the Cantor set 
and it further reduced asymptotically (as n — > 00) to a power law similar to the Gutenberg-Richter law for earthquake 
energies. In both the cases, the original distribution of the relevant variable (m 2 and S) was log-normal in which the 
normal factor became a very weak function of the variable in the asymptotically (t — ► 00 and n — > 00 respectively), 
thus rendering a power-law form to the distribution. Our derivations of the two power laws in the two vastly different 
models also indicate the universality of the exponents v = 7 = 1. In particular, the value of the Gutenberg-Richter 
exponent 7 — 1 in the fractal-overlap model is clearly independent of the dimension of the fractals used and therefore 
the result is of a general nature. In the context of this paper it may be mentioned that in a similar fashion Pietronero 
et al iTljj l found a common mode of origin for the laws of Benford and Zipf. 



G(S) ~ ~ for n 



00. 



(20) 



IV. CONCLUDING REMARKS 
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